1 

^  AD-A103  997 

OAKLAND 

UNIV  ROCHESTER  MICH 

F/G  12/1 

■I 

the  use 

OF  FOURIER  TRANSFORM  METHODS  FOR 

THE  EVALUATION  OF  COEF— ETC(U> 

■I 

AUG  81 

J  M  MCKINLEY »  P  P  SCHMIOT 

N00014-77-C-0293 

■1 

UNCLASSIFIED 

TR-6 

NL 

ntfc on  AD  A103997 


1 


The  Use  of 


Task  No.  NR  359-64^" 
Technical  ^eport^No .  6 


Wj  TP- 


Fourier  Transform  Methods  for  the  Evaluation  of  Coefficients 
in  a  Taylor  Series  to  Arbitrary  Order^ 


JO 

by)  J. 

'  D 


M . /McKinley 
P  ./Schmidt 


Prepared  for  Publication  in  the 
J.  Phys .  A  (Math  Phys.) 


Oakland  University 
Department  of  Physics 
Department  of  Chemistry 
Rochester ,  Michigan  48063 


Ji-  August  1981 


V 


1  0  1 3 SI 


> 


Reproduction  in 
any  purpose 


whole  or  in  part  is  permitted  for 
of  the  United  States  Government 


Approved  for  Public  Release;  Distribution  Unlimited 

81  9  10  010 


1  „  Rl  T  Y  CLASSIFICATION  OF  THIS  PAGE  'Wh  •"  Omtm  Entered) 


REPORT  DOCUMENTATION  PAGE 


1  REPORT  NUMBER 


4  T|  Tl  t  fond  5uMJfle) 


The  Use  of  Fourier  Transform  Methods  for 
the  Evaluation  of  Coefficients  in  a 
Taylor  Series  to  Arbitrary  Order 


Aj  THORfAj 

I  .  M.  McK in  lev 
P .  F.  Schmidt' 


HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


)  RECIPIENT'S  CATALOG  NUMBER 


i 


5  TYPE  OF  REPORT  ft  PERlOO  COVERED 

Techn i ca 1 
S/1/81  to  4/31/82 


6  PERFORMING  org.  report  number 


B  CONTRACT  OR  GRANT  NUMBER^ 

\ 0 0 0 1  I  -  77 -C- 0  239 


=  ERFORMING  ORGANIZATION  NAME  AND  AOORESS 

departments  of  Physics  and  Chemistry 
Oakland  University 
Rochester,  Michigan  48063 


’*  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


10  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


NR  359-648 


12  REPORT  DATE 


8/1/81 


13  NUMBER  OF  PAGES 

1  T 


AGENCY  NAME  ft  ADORES  Sf/f  different  from  Controlling  Office)  I  15  SECURITY  CLASS  (of  thte  report) 


o#  util  AERIFICATION  DOWNGRADING 
SCHEDULE 


’6  Distribution  statement  'of  this  Report) 

Approved  for  Public  Release 
1)  i  s  t  r  i  but  i  on  Un  1  im  i  ted 


*T  Distribution  ST  at  tMtN  i  tor  tne  aoatrect  entered  in  Woe*  JO,  tl  different  from  Report) 


«8  supplementary  notes 

Submitted  to  .1.  Phvs.  A  (Math  Phys.) 


19  KEY  WORDS  (Continue  on  reveree  aide  If  neceaeary  and  Identify  by  block  number) 

Taylor  Series,  Fourier  transforms 


20  ABf^b^ACT  'Continue  on  reveree  aide  If  neceeemry  and  Identify  by  block  number) 

We  present  a  new  formulation  of  the  Taylor  series  for  a 
class  oi  functions  which  can  he  expressed  as  the  product  of  a 
purely  angular  and  purely  radial  part,  viz., 


(:U>  =_vXl,( r)  f  lr) 


in  wliich  ( r )  is  the  spfftrrtcal  harmonic  function.  The 


20  cont inued 


\ 


formulation  is  based  upon  t  lie  use  of  the  lourier  transform  "I 
the  function  to  be  expanded.  The  general  expression  for  the 
Taylor  series  which  we  obtain  is  similar  in  form  to  the 
l.aplace  expansion  for  the  Coulomb  potential.  Thus,  Taylor 
series  can  be  developed  for  aribtrary  functions  along  lines 
which  are  similar  to  the  multipolar  expansions  in  the 
electromagnetic  theory.  Such  expansions  are  easily  adapted  to 
the  symmetry  of  a  collection  of  sources. 


'-1 


ffCuHlTY  CL  AMIAICATION  O *  fHI*  ►AGICWW  £>*•  t--rW 


The  Use  of  Fourier  Transform  Methods 


in  a  Taylor  Series  to 


for  the  Evaluation  of 
Arbitrary  Order 


Coefficient 


by 


J.  M.  McKinley 
Department  of  Physics 

Oakland  University 
Rochester,  Michigan  48063  USA 

and 

P.  P.  Schmidt 


Department  of  Chemistry 
Oakland  University 
Rochester,  Michigan  18063  USA 


- 1- 


Sunnnar  y 


IVc  present  a  new  formulation  of  the  Taylor  series  for  a  class 
of  functions  which  can  be  expressed  as  the  product  of  a  purely 
angular  and  purely  radial  part,  viz.. 


G(rJ  =  YAM(r)f(r) 


in  which  Y,  (r)  is  the  spherical  harmonic  function.  The  formulation 
is  based  upon  the  use  of  the  Fourier  transform  of  the  function  to 
be  expanded.  The  general  expression  for  the  Taylor  series  which  we 
obtain  is  similar  in  form  to  the  Laplace  expansion  for  the  Coulomb 
potential.  Thus,  Taylor  scries  can  be  developed  for  arbitrary 
functions  along  lines  which  arc  similar  to  the  multipolar  expansions 
in  the  electromagnetic  theory.  Such  expansions  are  easily  adapted 
to  the  symmetry  of  a  collection  of  sources. 


1.  Introduction 


Our  purpose  in  this  paper  is  to  present  a  new  and  powerful 
method  of  handling  the  expansion  of  a  class  of  functions  in  a  Tavlor 
series.  The  elements  of  this  class  arc  simply  those  functions  which 
can  lie  expressed  as  a  product  of  a  purely  radial  and  purely  angular 
part :  viz.  , 


C.fr)  =  VAii(r)  F  ( r  )  (1.1 

A 

in  which  (r)  is  the  spherical  harmonic  function.  We  proceed 
to  derive  the  general  term  of  the  Tavlor  series  to  any  arbitrary 
order.  We  are  able  to  show  that  given  the  Fourier  transform  of 
any  well-behaved  function  (or  generalized  function"),  it  is  possible 
to  obtain  the  specific  form  of  an  arbitrary  term  in  the  series. 

These  results  extend  and  generalize  an  initial  effort  which 
was  presented  elsewhere  by  us  (Schmidt,  Pons,  and  McKinley,  1980). 

In  particular,  by  more  laborious  means,  we  derived  specific  formulae 
for  the  first  and  second  order  terms  of  the  Taylor  scries.  The 
method  which  we  used  in  that  earlier  work  we  found  not  to  generalize 
easily  to  the  expression  of  a  general  term  for  the  Taylor  scries. 

The  results  of  this  paper  contain  exactly  the  results  reported 
previously  (Schmidt,  ct  al . ,  1980). 

We  find  several  reasons  for  undertaking  the  development  of  this 
new  expression  of  the  Taylor  series.  Foremost  is  our  need  to 
consider  the  displacement  of  a  particle  about  a  point  nf  equilibrium 
which  simultaneously  is  also  a  centre  of  symmetry  for  a  collection 
of  sources  (Schmidt,  et  al.,  1(18(1).  Thus,  we  need  to  consider 
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a  form  of  the  Taylor  series  expansion  of  an  arbitrary  function 
which  exploits  the  symmetries  of  the  system  in  the  same  manner 
as  is  customarily  done  for  the  electrostatic  potential  through 
the  use  of  the  Laplace  expansion  (c_f.,  Jackson,  1962).  These 
symmetries  are  most  effectively  handled  through  the  use  of  expansions 
which  depend  upon  the  spherical  harmonic  functions.  Hence,  it  is 
natural  to  consider  the  use  of  Fourier  transforms. 

Often,  and  for  molecular  physics  in  particular,  it  is  necessary 
to  consider  approximations  of  complicated  functions.  This  is 
especially  true  of  potential  energy  functions  which  are  derived 
from  molecular  integrals.  The  analysis  of  the  mechanics  of  molecular 
vibrations,  for  example,  commonly  depends  upon  approximations 
which  replace  the  exact  potential  energy  functions  with  harmonic- 
replicas.  It  is  universally  known,  of  course,  that  these  harmonic 
approximations  are  adequate  only  when  applied  to  the  analyses 
of  the  molecular  low-lying  vibrational  states. 

It  is  frequently  of  interest  to  examine  higher  order  terms  in 
order  to  assess  the  accuracy  of  the  harmonic  approximation.  In 
some  cases,  higher  order  terms  are  required  in  order  to  establish 
limits  of  stability  for  mechanical  systems.  For  a  relatively  simple, 
complete  function,  it  is  possible  to  extract  cartesian  harmonic  and 
higher  order  terms  as  individual  terms  in  a  Taylor  series.  Realistic 
potential  energy  functions,  however,  often  contain  fairly  complicated 
angular  dependencies.  These  angular  dependencies  reflect  the 
complicated  character  of  the  environment  which  surrounds  a  particular 
particle  (atom  or  molecule)  of  interest  (see,  for  example,  Uriels 
(1980)].  For  these  functions,  the  forming  of  the  Taylor  series 
ceases  to  he  a  simple  matter. 
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The  usual  representation  of  the  Taylor  series  presents  the 
determination  of  the  coefficients  as  a  sequence  of  differentiations. 
These  operations  are  applied  to  the  function  the  expansion  of  which 
is  required.  It  is  possible  in  principle  to  make  repeated  use  of 
Rose's  (1957)  formula  for  the  gradient  to  determine  these  coefficients 
in  terms  of  the  basis  vectors  of  the  spherical  tensors.  In  reality, 
however,  any  actual  attempt  to  apply  that  formula,  even  to  the 
simplest  of  functions,  quickly  becomes  unmanageable.  It  is  tedious 
even  to  determine  the  harmonic  terms  for  the  series. 

In  contrast,  we  are  able  to  obtain  the  coefficients  of  the  Taylor 
series  after  the  evaluation  of  only  a  single  form  of  radial  integral. 
The  strength  of  the  method,  therefore,  is  a  transparent  compactness 
and  flexibility  in  terms  of  ease  of  application.  This  strength  is 
illustrated  hv  several  examples  in  a  separate  paper  (McKinley  and 
Schmidt ,  198 _  b) . 

The  outline  of  this  paper  is  the  following.  In  section  2 
we  consider  the  single  centre  expansion  of  a  function  into  a  Taylor 
series.  A  particular  form  of  our  result  replaces  differentiations 
with  integrations.  For  complicated  functions,  which  possess  well- 
defined  Fourier  transforms,  these  integrals  prove  easier  to  consider 
and  evaluate  than  is  the  case  of  the  consideration  of  the  multivariate 
differentiation  of  the  function  directly. 

In  section  3,  based  on  the  use  of  the  results  of  section  2,  we 
present  a  derivation  of  Rose's  (1957)  formula  for  the  gradient 
operation.  A  consistency  between  Rose's  work  and  ours  is  established. 
In  addition,  we  use  the  differential  raising  and  lowering  operators 
for  the  spherical  Bessel  functions  to  re-establish  a  differential 
form  for  the  Taylor  series  coefficients.  This  differential  form 
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is  consistent  with  the  angular  dependencies  of  the  integral  form 
and  with  the  symmetries  of  the  spherical  harmonic  functions.  It 
is  clear  from  the  differential  form  obtained  that  in  many  instances 
the  integral  representation  is  indeed  the  simpler  way  to  approach 
the  determination  of  coefficients  in  a  Taylor  series. 

In  section  4  we  consider  the  Taylor  series  as  an  expansion 
about  two  centres.  It  is  possible  to  consider  the  general  develop¬ 
ment  of  any  well-behaved  function  as  an  expansion  about  two  or 
more  centres.  The  Carlson-Rushbrooke  (1950)  expansion  of  the  free- 
space  Green  function  is  such  an  example.  In  section  4  we  establish 
a  parallel  treatment  of  the  Taylor  series.  We  limit  our  considerations, 
however,  only  to  two  centres.  In  this  section  we  also  show  that  as 
the  two  centres  coalesce,  the  resulting  form  is  consistent  with  the 
single  centre  form  of  the  Taylor  series. 

Finally,  in  the  last  section  we  mention  the  limiting  forms 
for  scalar  functions.  Whereas  for  the  genera]  case,  for  a  general 
angular  dependence,  the  integral  form  of  the  Taylor  series  seems 
to  be  the  simpler  route  to  follow,  for  a  scalar  function,  we  find 
the  resulting  differential  form  is  equally  simple,  if  not  simpler,  to 
use.  Our  form  of  the  Taylor  series  in  its  differential  representation 
for  a  scalar  function,  however,  does  not  hear  a  close  resemblance 
to  the  cartesian  differential  form.  Oil”  form  preserves  the  angular 
dependencies  as  arguments  of  the  Legendre  polynomials.  Differentia 
tion  is  applied  only  to  the  radial  part  of  the  function. 
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2 .  The  lay  lor  Ser  i  cs  for  the  Expansion  about  one  Centre 


We  develop  the  Taylor  series  for  functions  which  are  separable 
into  purely  radial  and  purely  angular  parts,  viz.. 


G ( r )  =  F(r)YAu(r) 


(2.1) 


where  (r)  is  the  spherical  harmonic  function.  Still  more  general 
functions  could  be  expressed  as  superpositions  of  such  functions. 

The  function  G(r)  admits  a  Fourier  transform: 


G  ( r )  =  C  2  tt  ) 


-  j 


d  3 k  f  (k)YX)j  (k)exp(- ik*r) 


where  f ( k )  is  given  by 


f  Ik) 


4  tt  i 


dr  r2p(r)jA(kr) 
o 


and  j  A ( x )  is  the  spherical  Bessel  function  of  the  first  kind. 

The  function  at  the  displaced  point  r  +  c  is  given  by  the 
vectorial  Taylor  series  (Band,  1959) 


C(r+c)  =  l  X  f  c  •  V  ) n  G  ( r )  .  (2.4) 

n  =  0  ~  ' 


Consider  now  the  general  term  in  this  series  with  G(r)  given  by 
its  F . t  .  : 
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( c  -  V )  C  '.  (,  r ) 


1 


(271")  3n ! 


d  3  k 


f ( k ) V  ^  ( k ) ( - i c • k ) ncxp ( - i k • r ) 


J3k  knf(k)Y,  (k)  (c  •  k )  "exp  ( -  i  k  •  r )  . 
C2tt)  3n!  J  'Uj 


Here,  c*k  is  the  cosine  of  the  angle  between  the  unit  vectors  c 
The  nth  power  of  this  quantity  can  be  expanded  in  Legendre  poly¬ 
nomials  of  the  same  argument  (Morse  and  feshback,  1953) 


'N  ZN  /N 


The  coefficients  A  ^  arc  given  by 


A 


2£+l 


n£ 


r+l 


-1 


dx  x  P£(x) 


0,  5L > n  or  n-£  odd 


( 2  £  + 1 ) n !  (n - £  + 1 )  !  !  „  „ 

'  (n-£  +  i)~!  (n+£+l7!T’  *-n  and  n '  *  CVcn 


Wc  now  introduce  the  Rayleigh  expansion, 


exp  ( -  i  k  •  r )  =  [  ( -  i ) l'  (  2L+ 1 )  j  .  (kr )  1>  (k  •  r )  , 


and  the  addition  theorem  for  the  spherical  harmonic  functions, 


!Vk'rl  21, +  1  l -  VLMfk),>  LM(r)  ’ 

M 


2.5) 

a  n  d  k  . 


2.(0 


(2.7) 


(2.8) 


so  that  we  have 
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c\p  ( -  i  k  •  r )  --  4n  (  -  i  )  i  |  (  kr  (  Vj  M(k  |  Vj  M(  r 
L  ,  M 


(Ck)n  =  4,  l  (2^1)-1An!Ytm(k)Ytm(c). 
Z  ,m 


(2.11) 


With  these  substitutions,  Equation  (2.5)  becomes 


jn-(c-YV<Hr1  -  l  (-i)L*"(2t*n-',Vn(V|M(r)V,m(c) 

L  ,  M  ,Z  ,m 


*  *  * 


'  Jk  k-  Umkk,,f(k)jL(kr)YXy(k)YIM(k)Y5m(k) 
3  4  - 


(2.12) 


The  angular  integraL  is  the  well-known  integral  over  three  spherical 
harmonic  functions  (Rose,  1957) 


i  J^'\.,fk)YrM(k)Ytm(k)  = 


I  l  ‘  U  1I.UI  I, 

J  k  A u  EMV  ’  £m 
4  TI 


( 1 . P  0 0  |  AO)  (l.C.Mm  |  Ay1) 


(2.15) 


where  (I.CMm|Au)  is  the  Cl ebsch - Gordan  coefficient.  We  also 

introduce  a  quantity  defined  by  the  remaining  radial  integral: 


I,  (r)  =  — —  dk  kn+2f(k) j.  (kr) 
nL  (Ztt)  3  J  o  L 


(2.14) 


Altogether,  we  have  found  the  following  expression  for  the  general 
term  in  the  Taylor  series: 


(c  •  V )  (4  ( r )  =  ( 4 7T ) 


5/2  cn 


1  '-n,'*nA„,YrM(ovtm(c) 


L  ,M£  ,m 
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2L+ 1 


uwn  ( 2  x +1  ] 


1/2 


(L£00 | XO) (L£Mm| Xu) Tnl  (r)  . 


(2.15) 


The  most  frequently  used  cases  will  he  those  of  the  lowest 
order.  For  n=l,  £  can  only  be  1  and  A^  =  l.  We  have  then 


c-VG(r)  =  (4i»3/2c  l  (•i)L*1V]M(;)Ylm(o(j|l^1T 

L  ,  M ,  m  '■ 


1/2 


(L100 | X 0) (LlMm|Xu)  Ij,  (r) . 


(2.16) 


A  special  application  of  this  is  to  find  an  equilibrium  condition, 
which  can  be  expressed  as  c*VG(r)  =  0,  for  any  c,  or  as 


0  =  l  (-i)LYIM(^)(2L+l)1/2(L100|X0)(LlMm|Xu)I1I (r) .  (2.17) 

L,M  iL 

For  n=2,  l  can  be  either  0  or  2  with  A2q=1/3  and  A?2=2/3.  We  have 
then 


Lcc.vUcq)  -  c>fi(-i)**2YXii(hv 


0  0 


t  I  (•i)L*2Yll.(r)Y,mtc) 


L  ,M,m 


LM^  ;  2m 


2L  + 1 


5 ( 2  X  + 1 )  j 


* 1 2  L  ^  r  ) 


(C) I 2X (T) 

i  /2 

( L  2  0  0 | XO) (L2Mm| Xp) 

(2.18) 


These  lowest  order  terms  are  obviously  easy  to  evaluate,  but  even 
the  general  term  [Equation  (2.15)]  is  not  intrinsically  more 
difficult.  This  simplicity  traces  back  to  Equation  (2.11)  which  is 
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an  especially  simple  example  of  the  expansion  of  an  invariant  into 
irreducible  forms  (F ano  and  Racah,  1959). 

3.  The  Formula  for  the  Gradient 


In  this  section  we  verify  that  the  first  order  term,  eqn  (2.16), 
is  equivalent  to  the  usual  gradient  formula  (Rose,  1957). 

We  begin  by  evaluating  the  Glehsch-Gordan  coefficient  ( 1. 1 0  0 | >  0 ) 
in  expression  (2.16).  This  coefficient  restricts  the  parity  of  I. 
as  well  as  limits  its  range.  We  have 


1/2 

( A  -  1 ,100  |  AO)  = 


(A  +  l,  J00|A0) 


A  + 1 1 1 7  2  _ 

(tt+tJ  ' 


(3.1  ) 


all  other  values  of  L  yield  zero.  We  now  write 


c  •  VG  (r)  =  (-i)*^71^ 


3/2 


l 


2T+1 


l/: 


/3  M  ,m 


(A-l,lMm|An)YA_1M(r)Ylm(c) 


XIlfA-l<r>  +  [tXTtJ  (^MMm|Au)Yx  +  lMfr)Ylm(c)Ii>x  +  i(r) 

(3.2) 


We  handle  the  angular  and  radial  factors  separately.  The  standard 
definitions  of  the  vector  spherical  harmonic  functions  are  (rose, 
1957) 


I»,A*l,utr) 


J  IM(r)c 

M  ,m 


(3.3) 


where  the  quantities  f;  are  the  basis  vectors  for  the  spherical 
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tensors  (Rose,  1957).  The  scalar  product  of  (7.3)  with  the 
vector  c  is 


J  fA^Mn|Xp)Y  (r)c^ 
M  ,in 


=  /4ti7"3  l  (A±llMm|Au)YA  +  1  (^)Y  lc). 
M ,  m 


(3.4) 


On  the  other  hand,  for  the  radial  part  we  consider  the 
following.  In  the  definition  of  I^+^(r)  we  use  the  standard 
operators  for  lowering  and  raising  the  index  of  the  spherical 
Hcssel  functions  (Morse  and  Feshback,  1953,  and  Infeld  and  Hull, 
1951).  Thus, 


dk  k  3  fOO.i  x_ ,  (krl 

) 

o 

dk  k3  f (k) 


J  o 


(ill  +  d 

(~r~  HT 


TF-Vkr)  *  (TTTTT  h(kr) 


dk  k2  f(k) j A(kr) . 


(.3.5) 


Similarly, 


(2^)3In  +  l(r)  =  j  dk  k3f(k).iA  +  ]  (kr) 


A  _  d 
r  efr 


dk  k2  f (k) j A (kr) . 


(3.(0 


At  this  point,  we  identify  the  inverse  of  eqn  (2.5) 


Dr) 


_  (- i ) 


2tt  2 


dk  k2f(k)jA(kr) 


(  3 .  "  ) 
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Upon  subst i tut  ion  of  all  of  these  expressions  into  eqn  (3.2), 
we  find  exactly  the  standard  formula  for  the  gradient  operation 
{Rose,  1957) : 


This  new  form  in  terms  of  Fourier  transforms  could  have  been 
anticipated,  because  the  radial  operations  in  eqn  (3.8)  are  uniquely 
those  associated  with  the  spherical  Bessel  functions. 

Fqn  (3.8)  completes  the  connection  with  Rose's  formula  for 
the  gradient.  We  now  extend  the  analysis  given  above  to  the 
consideration  of  an  arbitrary  term  in  the  Taylor  series.  Thus,  we 
consider  the  transformation  of  an  arbitrary  term  from  an  integral 
to  a  differential  form  in  the  radial  component  only. 

For  a  given  function  G(r),  A  is  specified.  For  a  given  value 
of  the  order  n,  the  values  of  SL  also  arc  specified.  The  Clebsch- 
Gordan  coefficients  fix  the  parity  and  range  of  tiie  L-values. 

Thus,  it  is  possible  to  see  that 

b  =  A  +  q  (3.9) 


with  the  values  of  q  fixed  and  limited  by  the  conditions  mentioned. 
The  summation  over  I,  in  the  expression  for  an  arbitrary  term  in  the 
Taylor  series  is  rewritten  as 
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(i*-V)nr.(r)  =  1471) 


3/2  c 


1  (-i)"** — 


11! 


n ! 


ne 


[  (2c+i)  (2d+D  ] 


\  .1/2  e.m 


v0„.(c) 


X  I  I 

q  =  0  M 


( -  i ) q  /2  (A+q-)  +T ( A+q  £00  |  X0)  ( A+q  £Mm  \  Xu)  YA  +qM  ( r  )  I„  ^  +q  f  r ) 

A 

+  (-i)-q/2(A  -  q  )'+T  ( A  - q  £0  0  |  AO)  (  A -q  £Mm  |  Ay )  Yx  _qM  ( r )  In  ^  A  _q  (r ) 


(3.10) 


In  fact,  the  summations  run  over  q=0,2,4,6...  or  q=l,3,S,... 
depending  upon  the  values  of  n  and  A.  The  factor  together  uith 

the  Clehsch-Gordan  coefficient  automatically  sort  out  the  terms 
which  survive  in  the  summation.  Hence,  there  is  no  need  to  display 
any  specific  form  of  the  sorting  process. 

We  now  use  the  following  lowering  and  raising  operators 


Jx+q<kr> 


=  (-D 


.A  +  q 


rdr 


q  -  A  r  i.  ~  A  . 
r  (k  j 


(kr)) 


(3.11) 


and 


JA-q^  =  r 


q  -  A  -  1 


d 

Fd7 


■A  +  1{k'qiA(kr)l 


(3.12) 


to  write 


<2’>’In,AVr)  '  (-|)t,r>‘q 


f  d  lq  -  A 

rilrj  r  J 


dk  kn"q+2f(k)j A(kr)  (3.13 


and 


<2’>’In,J-<l(r’  ‘  r" 


-A  -  1  (  d  VLa  +  1 
rdr 


dk  kn'q+2f(k) jx(kr)  (3.14) 
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For  a  given  value  of  q  in  T  .  (r)  we  have  ad i listed  the 

n  ,  A  7  q 

index  on  the  Bessel  function  to  a  value  of  X.  It  remains  therefore 
only  to  lower  the  magnitude  of  the  exponent  of  k  in  the  integrand 
from  k11  q  +  “  to  k“.  This  we  do  as  follows.  First,  note  that  from 
the  various  conditions  on  the  parities  of  the  indices,  n-q  will 
always  be  even.  Thus,  we  can  reduce  n-q  to  a  zero  value  in  steps  of 
2.  From  the  eigenvalue  equation  which  the  Ressel  function  i ^ 
satisfies,  we  write 


x2.i  A  (x) 


A(A+1)  '  ^rfx2^)]jx(x')- 


(3.15) 


The  right  hand  side  of  this  equation  defines  an  operator  which  lowers 
the  exponent  on  k  (x=kr)  by  2.  The  use  of  this  operator  (n-q)/2 
times  yields  the  desired  result: 


dk  kn"q+2f (k) j. (kr)  = 


n-u 


r 


dk  k  f  ( k )  i  ^  ( k  r  1 


2  tt  2  i 


.  X 


n-q 

d  r  _2  d  J  2 


A(A  +  1)  '  ^TF(r2dT)  rq_nF  (r) . 


(3.1b) 


The  final  expression  for  the  arbitrary  term  in  the  Taylor 
series  is 


n 


i  (cV)n(l(r)  =  /4tt  l  (-i)n - 

n-  ~  ~  ~  n>  ?.,m  [  ( 2  X  + 1  (22  +  1)1 


T77  YPm(cl 


■  l  l  iq 

q=  0  M 


/2(X+q)+l  ( X  +  q  5, 0  0  j  XO)  (X  +  q2Mm|  X  P )  Y  ^  M  ( r )  r 


X+(l 


d 

rd  r 


l<  -A 

r 


♦  i 


2(X-q)  +  l(i-qfOO|XO)  (X-q2,MmlXp)YA_qM(r)rq'A'1  f^]  'r>  +  1 


TVAT1J  '  ^TFir2^FJj 


(3.17) 


t.  The  Expansion  ol  a  Function  in  a  Tavlor  Series  about  two  Centres 


tor  some  applications  it  is  useful  to  consider  the  function  of 
ci|ii  (2.1)  to  depend  upon  two  points  ri  and  r2  ,  cither  of  which 
in  ay  be  displaced  independently.  We  set 


r  =  r,  -  r ; 


(1.1) 


and  redefine 


b(r,  ,  r  2  )  =  Y^l  r)  l;(r)  . 


(■»•••> 


Ibis  function  still  possesses  a  Courier  transform  as  expressed  by 
eqn  (2.2),  (2.3)  and  (3.7).  Before  we  evaluate  this  function  at 
displaced  points,  however,  we  present  an  alternative  Fourier  transform 
which  emphasizes  the  two-centre  character. 

Into  the  expression  for  the  Fourier  transform 


C(rl,r2)  =  (2n)  '  d3k  f  (k) Y ,  ( k) exp ( - i k • r i ) exp ( ik* r 2 ) 


(4.3) 


we  substitute  tmo  Rayleigh  expansions  of  the  form  (2.10).  The 
angular  integrations  over  the  three  spherical  harmonic  functions 
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G(r i , r2 ) 


fo 


dkk' 


4  71  Lj  ,  Mj  ,  L2  ,M2 


:jLl (kr1)YLiMi Cri)YLlMl (k) Jl2 (krz)YL2M2 (r2)YL2M2 


3  ~H> 

it  '  LiMiL2M2 
1/2 


I  iL2'L,YLlM/ri)Wr2) 


X  C.lL.i.^l.LpLWj  (LiL200|A0)(L,L2MiM2|Au) 


r 

dk  k2  <  (k)  j  Li  (kri  )  jL?  (kr2 ) .  ( 4 . 4  J 


In  order  for  eqn  (4.4)  to  be  consistent  with  eqn  (3.7)  and  (4.2) 
for  arbitrary  functions,  we  require  the  following  double  addition 
formula  for  spherical  harmonic  and  Bessel  functions: 


Y.  (r)j.(kr)  =  (4tt)  1/2  £  i  A  +  L2  'L  1  (L  i 1. 2  00  |  AO)  (L  !  L2M  ,  M2  |  Ap) 

Am  A  L iMi L2M2 


(2L i + 1 ) (2L2+1)1 

zrn 


1/2 


YI.1M1^ri)jL1(kri)YL2M2(r2)jL2Ckr2)-  (4-M 


A  similar  form  for  real  spherical  harmonic  functions  has  been 
derived  by  Johnson  (1973).  It  is  possible  to  show  easily  that  the 
right  hand  side  of  (4.5)  reduces  to  the  form  of  the  left  hand  side 
when  r2=0  and  r=ri.  Additionally,  when  A=y=0,  the  right  hand  side 

/N 

of  (4.5)  can  be  shown  to  reduce  to  Y 0 o (r) j o (kr)  with  the  use  o(  the 
appropriate  Clcbsch-t’iordan  coefficients  and  the  addition  theorem 
for  the  spherical  harmonic  functions.  Thus,  the  two-centre  expansion 
is  consistent  with  the  more  familiar  single-centre  forms. 

We  proceed  to  evaluate  the  two-centre  function  at  the  displaced 
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points  i']+a  and  r2+b.  We  have  as  the  extended  vectorial  Taylor 
series,  the  expansion 


G ( r  j  +  a , r 2  +  b )  =  £  -^-(a-Vi  +  b • V2 ) nG (r j  , r 2 )  . 

~  ~  '  ~  n=0  n-  ~  ~  -  ~  ~  ' 


( 4  .  (> ) 


The  general  term  in  this  series  is  to  be  evaluated  with  G(ri,r2) 
expressed  by  its  Fourier  transform  [eqn  (. 4 . 3 J  ]  .  We  have 


jYy(.a*Vi  +  b*V2)  G  (r  i  ,  r  2  ) 


=  (2tt)'3| 


d3k  f (k) VAu (k) 


x  ppy[-ia*k  +  i b  •  k  J 11  exp  ( -  ik  •  r  i )  exp  ( ik  •  r  2 )  . 


(4.7) 


In  this  expression  we  use  the  binomial  expansion  and  the  expansion 
in  irreducible  forms  (2.11)  to  replace  the  quantity  in  the  brackets: 


[a-k  -  b  -  k  ] 11  =  E("Ha-k)n'q(-b-k)q 


y  r 11 1 1. 11  -i n ~ q r  hiq.  C411)2 ^ 

o  l  qk  (  b]  (7TT+T7 (2I7TTT 

q  ,  £ ,  ,  S,  2  m  i  ,  m  2  1  *•  v 


x^n  -  q ,  £  i  Aq  ,  £  2  Y  £  i  m  i  £2m2  £2in,  ' 


(4.8) 


The  coupling  rule  for  the  spherical  harmonic  functions  is  now  used 
to  reduce  the  number  of  these  functions  in  the  product  (Rose,  1957): 


is  ''•ft  a 

Y£1m1  (k)Y£2m2(kJ  =  TO  *+11  J 


x  (£1  £?ni,m2  I  £m)Y£m(k)  . 


(1.9) 
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The  substitution  of  two  Rayleigh  expansions,  and  the  use  again  of  the 
coupling  rule  together  with  eqn  (4.8)  and  (4.9)  leads  to 


ip-[a*Vi  +  b*V2]nG(r1,r2)  =  1  I 

q,£i,H2,mi,in2,Li,L2,Mi,M2,L,M 

"(-1)n*L,-L2(;)an-<'(-b)X.q>lliAq  o)Y,2raia)YLlHl(n) 

1/2 

xYL2M2,;r2^  T2  £  i  + 1 )  {  21 2 +11  ('2  F+1TT2 LTiy)  Uifi-aOO  |  £0)  (£i£2»i1m2  |  Cm) 


x (L i L2 00 | L0)  (Li L2MiM2 | LM) j  dk  kn+2 j  ^ (kr j ) j j ^ (kr2) 


4n 


dV\p^YLM^Y£m^ 


(4.10) 


The  angular  integral  is  evaluated  as  before  (cf.,  eqn  (2.13)),  and 
we  introduce  another  quantity  defined  by  the  remaining  radial  integral: 


f 00  ^ 

J  ,  ,  (r i »r2)  =  — - —  dk  kn+2f(k)j,  (krx)j  (kr2). 

nLiL2  -  ~  (2-n)  3  J0  Li  L2 


(4.11) 


Thus,  we  obtain  the  following  expression  for  the  general  term  in  the 
Taylor  expansion  of  a  two-centre  function: 


3/2 

jr  i  t  ^  *  V  i  +  b.V2]nG(r!  ,r2)  =  -  I 

n !  q,£i,£2,mi,m2,Lj,L2,Mi,M2,L,M 


1/2 

U.^OOUO)  (L,Lj00|L0) 


(fl.OOj  AO)  (£1£2m1m2|£m)  (LiL2MiM2  j  LM)  (£LmM|  Ap).JnJ  L  (rltr2).  (4.12 
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A1  tomato  expressions  lor  the  expansion  about  two  centres  are 
possible.  We  present  the  formula  for  the  general  term  of  the  Taylor 
series  for  one  such  modification.  In  particular,  when  the  vectorial 
difference  ri  -  r 2  =  R,  a  vector  drawn  between  two  points  of  reference, 
is  a  constant  quantity,  then  the  expansion  can  be  considered  in  terms 
of  displacements  about  the  two  ends  of  the  vector  R.  The  formula 
for  the  general  term  is 


L  [a-Vi  +  b •  v 2 ] ni; (R)  =  i^r-2  l  (-i) 

q,J.  j  ,  ,m.  ,m2  ,£,1,,M 


n  +  L 


xln)an‘q(-b)qA  .  A  . 
V  n-q,£l  q,H2 


21,+ 1 

[{2i  rn )  ( 2 1 2  ni  r  2  a+tt 


1/2 


(l!,«  200|  iOHUl.OOl  AO)  lt,*.2m,m2  |£m)  Itl.mMl  Ap)V  la)VJ<imj  (b) 


xVLmIR>WR> 


(1.151 


in  which  1  .  (R)  is  defined  by  cqn  (2.14).  The  analysis  which  leads 
to  eqn  (4.15)  parallels  the  proceeding  analyses. 


5.  Discussion 


Many  specific  examples  of  functions  which  can  be  expanded 
usefully  in  Taylor  scries  easily  spring  to  mind.  Several  examples 
which  can  be  derived  from  the  Yukawa  potential  arc  considered  in  a 
companion  paper  (McKinley  and  Schmidt,  198_a). 

Wc  conclude  this  paper  by  showing  the  limiting  forms  for  the 
Taylor  scries  when  the  function  to  be  expanded  is  scalar. 


1 


The  scalar  function  H(.r)  is  expressed  as 


G(.r)  =  Y  o  o  ( r)  [/4?F(r)] 


=  FCr) . 


The  Fourier  transform  of  this  function  is  Y00(k)f(k)  with  f(k) 
given  hv 


(5.1) 


f(k)  =  14 n) 


f  tx> 

3/-[  dr  r 


F(r) j o (kr) . 


(5.2) 


In  terms  of  this  form  of  scalar  function,  the  general  term  for  the 
Taylor  series  is 


ffl  ■  (4-)3/2(cn/n!)  I  (•i)n*l,A  jYj  (c)VLM(r) 

l ,  m ,  L ,  M 

*/ (2L  +  lj  (2£+TJ (hS-00  I  00)  (L£.Mm |  00)  I  j  (r) 


,  .  ',3/2,  n,  , v  r  -\n+£A  ^  ^ 

*  (4 n)  (c  /n!)  I  C-i)  A  . - 


-  >nl^ 


=  /4tt  (cn/n!  )X(“i)n+  £An£P  k.  (r-c)  lnP  (r)  . 


(5.3) 


This  quantity  is  evidently  a  generalization  of  the  expansion  of  the 
Coulomb  potential  in  a  Legendre  polnomial  series: 


r^d'1  =  lVr*c) 


The  general  differential  form  also  assumes  a  simpler  appearance 
in  the  scalar  limit.  We  need  only  consider  the  ln^(.r)  quantity  in 
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eqn  (5.3).  !  rom  the  definition  of  1  (r)  we  have 


dk  kn  +  “ f (k j  j  p  fkr ) 


=  4,.5/2(-l)<n+t>/2r4 


d 

rTTr 


l  (d/dr]n"4rl;(r) 


(5.5) 


following  Todd,  Kay  and  Silverstone  (1970),  we  can  show  that 


t  r  d 

rJr 


r'1  (d/dr)n‘4rF(r) 


y  .  f-l)l|(^q)  ! 
q  =  0  f  ii-qT!  (2q)  !  ! 


n  -  q  ,  d 
r"  dr 


x(d/drjn'q_1l;lr) 


(  5.0) 


Thus  , 


n  v 


1  i  j  ( n  ■*■£)/  2  f  (-l)q(£  +  q)! 

Ai  '  Jo  IKiTWTI 


-q 


n  -  q  d 

r  J7 


(d/dr)n‘q'lF(r). 


15.71 


for  example,  the  second  order  term  in  the  Taylor  series  assumes 
a  s i mp 1 e  form :  vie., 


i  2  di  +  d2r 
3lr  dr  dr 2 


»P i ( r • c ) 
.•>  • 


1  dl 
r  d  r 


d 2 1 
d  r ?  > 


(  5  .  S  ) 


This  form  is  particularly  useful  loi  a  number  of  applications  which 
involve  scalar  functions. 
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